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We present spatial and dynamic information on the s = 1/2 distorted kagome antiferromagnet 
volborthite, Cu3V207(OD)2 • 2D2O, obtained by polarized and inelastic neutron scattering. The 
instantaneous structure factor, S{Q), is dominated by nearest neighbor pair correlations, with short 
range order at wave vectors Qi = 0.65(3) and Q2 = 1.15(5) emerging below 5 K. The 

excitation spectrum, S{Q,(jj), reveals two steep branches dispersing from Qi and Q2, and a fiat 
mode at cj/ = 5.0(2) meV. The results allow us to identify the cross-over at T* ^ 1 K in ^^V NMR 
and specific heat measurements as the build-up of correlations at Qi. We compare our data to 
theoretical models proposed for volborthite, and demonstrate that the excitation spectrum can be 
explained by spin-wave-like excitations with anisotropic exchange parameters, as also suggested by 
recent local density calculations. 

PACS numbers: 75.10.Jm, 75.10.Kt, 78.70.Nx 



The quantum kagome Heisenberg antiferromagnet 
(QKHAF) is among the most coveted targets in the 
quest for experimental realizations of quantum spin liq- 
uid ground states. The two most prominent physical real- 
izations of the QKHAF studied to date are the naturally 
occuring minerals herbertsmithite, Cu3Zn(OH)6Cl2 [1| 
and volborthite, Cu3V207(OH)2 • 2H2O [2]. The physics 
of herbertsmithite is arguably influenced by depletion of 
the kagome lattice caused by antisite mixing, possibly re- 
sulting in a valence bond glass state [l,!^- In volborthite, 
the kagome planes are slightly distorted, but the lattice 
coverage is essentially complete. Since both systems de- 
viate from the pure QKHAF model, focus has shifted to 
the intriguing question of which states arise when the 
QKHAF is perturbed. Remarkably, a multitude of dif- 
ferent states have been proposed theoretically, depend- 
ing on the nature of the perturbation. These range from 
ordered states to the aforementioned valence bond 
glass. Experimentally, the ground state of volborthite re- 
mains enigmatic, despite over 9 years of intensive study 
Here we present a neutron scattering investi- 
gation of volborthite, employing both xyz polarised and 
inelastic time of flight techniques. 

Volborthite contains distorted kagome planes of edge 
sharing Cu^+ octahedra, well separated (~ 7.2 A) by 
pyrovanadate columns (V2O7), Fig. [I{a). There are 
two crystallographically distinct Cu^^ ions. Cu(2) form 
chains along the b directions whileCu(l) populates the in- 
terchain sites. Locally, Cu(l) and Cu(2) reside in tetrag- 
onally and axially distorted octahedra, respectively, sug- 
gesting the <i3z2_r2 orbital being singly occupied on the 
Cu(l) site, with d^2_y2 the magnetically active orbital for 



Cu(2). As a consequence, there are two different nearest 
neighbour exchange pathways, Ji (rcu(2)-Cu(2) = 2.93 A) 
and J[ (rcu(i)-Cu(2) = 3.03 A). Ji links Cu(2) ions along 
the b axis. J[ connects Cu(l) and Cu(2) ions. Further- 
more, the edge sharing of the Cu(2) octahedra along the 
b direction imply the possibility of a strong next-nearest 
neighbour exchange, J2. 

Despite a large estimated average coupling Javg = 
(2J{ + ^i)/3 = 84 K, magnetic susceptibility and spe- 
cific heat measurements on volborthite show no signs of 
long range order down to 1.8 K [2|. At yet lower temper- 
atures, muon spin rotation (/iSR) and ^^V NMR loHlij 
indicate slowing down of fluctuations at T* = 1 K, but 
with dynamics persisting to 20 mK. The state below T* 
was interpreted as either incommensurate or short range 
correlated. Consistent with this, low temperature spe- 
cific heat studies indicate a high density of low energy 
modes below T* 

These results, however, have provided little insight to 
the nature of the magnetic correlations and excitations. 
Neutron scattering is an ideal probe for investigating 
these aspects, but has thus far not been employed, in 
main due to the large background generated by spin in- 
coherent scattering from ^H and ^^V. We minimised this 
problem by replacing most of the ^H (~ 98%) by ^D, 
which was achieved by performing the synthesis [2| us- 
ing deuterated reagents in an atmosphere of N2 , followed 
by annealing the product several times in D2O at 95° C. 
Phase purity of the resulting sample was verified by pow- 
der XRD. The concentration of paramagnetic defects, 
proportional to the magnitude of the Curie tail Q , was 
estimated to be to be < 1%. 
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FIG. 1: Structure of volborthite. (a) Kagome planes of Cu(l) 
and Cu(2) octahedra (respectively blue and green) are sep- 
arated by V2O7 columns (red), (b) Local environments of 
Cu^^ in a kagome plane. Solid and dashed black lines indi- 
cate the two nearest neighbour exchanges, Ji and J[. The 
next-nearest neighbour exchange J2 along the 6-direction is 
shown in orange. 
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FIG. 2: (dcr/dQ)mag for volborthite at T = 200, 15, 10,5 K. 
(a) (dcr/dQ)mag mimics /(Q)^ (dashed black line) for Cu^^. 
(b,c) Spectral weight shifts into a broad feature around 1.1 — 
1.4 Green line is fit to Eq. (1). (d) Fit to (dcr/dQ)mag at 

5 K as described in the text (red). The green line indicates 
the fit at 10 K, showing the shift of spectral weight into the 
two sharp features at Qi and Q2' 



Polarised diffuse neutron scattering was performed on 
D7 at ILL using 35.7 g of sample and incident en- 
ergy Ei = 8.5 meV. Three orthogonal neutron polari- 
sations (x, and z) and their corresponding spin flip 
and non-spin flip cross sections were analysed, which 
allowed for isolation of the magnetic scattering cross 
section, {da/dQ)mag As no energy analysis was 

used, the observed scattering was effectively integrated 
up to = 8.2 meV, thus approximating the instanta- 
neous structure factor, S{Q). Spectra were measured at 
T = 200 K, 15 K, 10 K, and 5 K (Fig. [2j). 

At 200 K~ 2 J, the Q-dependence of (dcr/dr^)mag ap- 
proximately follows the Cu^+ form factor, |/((5)P, as an- 
ticipated for a paramagnet. Reducing the temperature 
to 15 K, broad diffuse scattering develops around Q = 
1.1 — 1.4 A~^. The Q-dependence is consistent with a 
buildup of nearest neighbour pair correlations, described 
by the powder averaged structure factor {d(j/dQ)mag- 

(1) 

where the second term in the parenthesis reflects the aver- 
age correlation (So • Si) = —0.25(5) between a unit spin 
and its Zi nearest neighbours at a distance rcu— Cu 

3 A. 

The total scattering was found to be 0.99(8)/i^ per Cu, 
corresponding to 33% of the full S{S + l).The fact that 
correlations are weak and confined to only nearest neigh- 
bours even at T/ Javg ~ 0.2 are both indicators of strong 



frustration in volborthite. 

As T is further decreased to first 10 K and then 5 K, 
the broad diffuse scattering persists, but some (14%) of 
the spectral weight is shifted into two sharper (though 
not resolution limited) peaks at Qi = 0.65(3) A~^ and 
Q2 = 1.15(5) A~^. The corresponding correlation length, 
^ = 24(8) A- 8rc u_Cu7 was extracted by fitting the dif- 
fuse scattering to (1) and the two sharper features to 
Lorentzians (Fig. [2](d)). The dynamical structure fac- 
tor, S{Q^uj), was investigated by inelastic neutron scat- 
tering on IN4 at ILL. Spectra were collected using in- 
cident energy Ei = 17.2 meV, giving elastic Q range 
0.65 A~^ — 4.95 A~^. Further experiments were carried 
out on MARI at ISIS {Ei = 15 meV, 0.45 A-^ -4.95 A"^ 
at a; = meV). Data are summarised in Fig. [3l 

At 50 K~ J/2, S{Q^uj) shows only a broad response 
centered at Q = 1.1 A~^ and extending to 5 meV, con- 
sistent with fluctuations in a short range correlated sys- 
tem. The Q-dependence extracted by integrating over 
the range 1 — 6 meV is indicated in the top panel of Fig. 
[3Kb). Its form is similar to S{Q) at 10 K and 15 K, and 
can also be fitted by Eq. 1 using r ~ 3.5(2) A~ ^cu-Cu- 
Acoustic phonons are observed dispersing from nuclear 
Bragg positions at Q > 2 A~^ and intense phonon scat- 
tering is found above 7 meV, making extraction of the 
magnetic signal at these energies difficult. Cooling to 
5 K, the low Q, low uj intensity has largely moved into 
two features: an intense broad flat band, centered at 
Uf = 5.0(2) meV, and a nearly vertical bar of scattering 
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FIG. 3: (a) S{Q,uj) from IN4 at T = 0.05, 5, 50 K, and MARI at 1.7 K. The integration ranges used for the cuts in panels 
(b) and (d) are shown as black rectangles, (b) Q-dependence of intensity over energy ranges and 3-4 meV for 0.05-5 K (fit by 
two resolution convoluted Lorentzians) and 1-6 meV for 50 K (fit by Eq. 1). (c) Temperature variation of intensities of the Qi 
(open circles) and Q2 (closed circles) modes. T* = 1 K is indicated by a red vertical line, (d) Q-dependence of the fiat mode 
integrated over 4.5-5.5 meV. S{Q) for a singlet-triplet excitation is shown in red. 



at Q = 1.08(2) A~^, which coincides with the Q2 peak in 
S{Q). Both of these features sharpen as temperature is 
reduced towards 0.05 K, with a second bar of scattering 
at Q = 0.68(4) A"^ - Qi growing below 1.7 K. Q-cuts 
through the Qi and Q2 modes are shown in the lower 4 
panels of Fig. ^h). While the peaks narrow somewhat 
with decreasing T, a more dramatic change is observed 
in their respective intensities, I{Qi) and I{Q2)- On cool- 
ing, I{Q2) remains constant, while I{Qi) increases to a 
final ratio I{Qi)/I{Q2) = 1.6 at 0.05 K. The buildup of 
dynamical correlations at Qi position thus coincides with 
the transition at T* observed in NMR and /iSR. 

The lineshape and ampliude of the fiat mode, on the 
other hand, show little temperature dependence, with 
only slight narrowing to become resolution limited be- 
tween 5 K and 1.7 K. Such narrow fiat modes are often 
associated with two level excitations, eg. between a sin- 
glet and triplet. Indeed, such excitations are expected 
in the site diluted QKHAF 0, [Tej. While the predom- 
inantly nearest neighbour correlations observed in S{Q) 
are consistent with such a state, the Q-dependence of the 
fiat mode does not match the sine let-triplet S{Q) [1, 
[Fig. 3(d)]. Another possible explanation for the fiat 
band is thus that it is associated with the short range 
order observed at the Qi and Q2 positions. This leads 
to two scenarios: i.) That the ground state possesses a 
degeneracy which results in a fiat band, as is the case in 
the ground state of the classical kagome systems [l7| or 
ii.) That the fiat band arises as a consequence of powder 
averaging at the zone boundary of a spin wave dispersion. 

In the pure QKHAF, no short range order is expected, 



contrary to observations. To understand the low T short 
range ordered state, it is therefore necessary to look the- 
oretically further afield. From structural considerations 
outlined in the introduction, the Hamiltonian of volbor- 
thite can include three exchange terms (Ji, J(, and 
J2) as well as a Dzyaloshinskii-Moriya (DM) interac- 
tion, for which the dominant component is considered 
to be Dz. The subsets of this model which have been 
treated semi-classically or quantum mechanically are: 
Ji = J( > 0, J2 = 0, Dz 7^ 0, the isotropic kagome DM 
model (IKDM) 0,[i,[l8| and Ji 7^ J{ > 0^2 =0,Dz = 0, 



the anisotropic kagome model (AK) In addition, 

recent L(S)DA+U calculations have suggested a model 
where -J[/Ji - 1.2 - 2, -J2/J1 = 1.1 - 1.6, and 
J( = 8.6 meV. This model, which we call the coupled 
chain model (CC), has only been treated classically so 
far [l9| . The 5 possible ordered states which emerge from 
these models are: q = order (IKDM for L)^ > 0.1, AK 
for 1/2 < J[/Ji < 1), q = a/3 X (AK for 1 < Ji/J[ < 
1.3), chirality stripe (AK for Ji/J[ > 1.3), ferrimagnetic 
(AK for J(/Ji < 1/2, CC for J2 < | Ji|/4 + J(/8), and 
spiral (CC for J2 > | Ji |/4 + J(/8). The position of Bragg 
peaks for these structures are compared with our exper- 
imental Qi and Q2 in figure 4(a). None of the proposed 
structures yield strong scattering at Qi. 

Therefore, instead of attempting to describe the ob- 
served S'(Q,cj) in terms of the microscopic model above, 
we take the phenomenological approach of constructing 
a generic dispersion emanating from antiferromagnetic 
zone centers at high symmetry positions on the circles 
defined by Qi and Q2 (Fig. 4a). One example of such a 
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dispersion among the scenarios considered assumes zone 
centers close to the (10) and (01) positions, yielding the 
dispersion: 



w(q) = ^(2J| + 2Jb^)2 - [2J|cos(fea) + 2J'^cos{qyb)]^ 

(2) 

where ^ are effective exchanges, giving the amplitudes 
of the dispersion along a and b. Then, 



S{ci,uj) = \F{r)\ 



J 2 — cos{qxa) — cos{qyb) 

^(q) 



(3) 



where \F{r)\'^ is the structure factor at the chosen po- 
sitions, 1/uj describes antiferromangetic spin wave in- 
tensity, and the numerator is a geometric term yielding 
zero intensity at ferromagnetic zone centers. To yield 
a smooth continuous function, 6'(q, cj) was interpolated 
between adjacent Brillouin zones. Finally, the spectrum 
was powder averaged and convoluted with the experi- 
mental resolution. The result of this procedure using 
J| = 5.1 meV and = 15.3 meV closely resembles 
the experimental data (Fig. 4(b)). One important con- 
clusion from this analysis is that the steepness of the 
Q2 mode implies that the 5 meV band cannot be the 
the global zone boundary energy. Instead the flat band 
is found to be a saddle point, which requires sizeable 
anisotropy between J| and J^. This result is consistent 
with LS(D)A+U calculations, which also suggest a sig- 
nificant anisotropy in exchange along a and b. 

In summary, we have reported polarised and inelastic 
neutron scattering results on the quasi-kagome 5 = 1/2 
antiferromagnet volborthite. These reveal three impor- 
tant features: i) Buildup of nearest neighbor pair correla- 
tions from 50 K to 10 K; ii) Short range order indicated 
by peaks at Qi = 0.65(3) A-^ and Q2 = 1.15(5) A-^ 
in the diffuse and inelastic scattering below 5 K; and iii) 
The excitation spectrum displays dispersive modes em- 
anating from both Qi and Q2 along with a fiat mode 
at cj/ = 5 meV. The inelastic intensity at Qi becomes 
dominant below 1.7 K, identifying the nature of the low- 
T state reported from ^^V NMR and Cp. Most mod- 
els proposed for volborthite however inaccurately predict 
no scattering at Qi, and could therefore be ruled out. 
We provide an empirical dispersion model which repro- 
duce the experimental observations for an anisotropic set 
of exchanges which are in rough agreement with recent 
LS(D)A+U calculations. A quantum treatment of this 
model could well yield a correct description of both the 
correlations and excitation spectrum that we have un- 
covered experimentally. 

We thank Z. Hiroi, B. Fak, P. Mendels, F. Bert, and 
O. Janson for stimulating discussions, G. Ehlers and J. 
Gardner for assistance during initial experiments, and 
SNF for support. 
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FIG. 4: (a) Reciprocal space of volborthite: the structural 
unit cell is indicated by the rectangle, and the extended Bril- 
louin zone of the kagome lattice by the dotted hexagon. Sym- 
bols represent Bragg peaks of the orders listed in the text, 
(b) The experimental S{Q,uj) measured at 0.05 K compared 
with the powder averaged S{Q,(jj) derived from our empirical 
spin wave model. The dashed line in the left panel indicates 
the (Q,c<;) window of the experiment. 
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